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Theoretical progress on \Vus\ on lattice 



1. Introduction 

Accurate knowledge of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element |Vi„ | is im- 
portant because it gives the basic parameter A in the Wolfenstein parametrization of the CKM 
matrix and is relevant to a stringent test of CKM unitarity + \ Vus \ + \ Vub\^ = 1. The K12 and 
Ki3 decays provide two precise determinations of |V„.v|, where their dominant uncertainty originates 
from theoretical evaluations of hadronic matrix elements, namely fx/ fn and Kij, form factors. 

Lattice QCD can provide a non-perturbative estimate of these matrix elements from first prin- 
ciples. Due to the limitation of the computational resources, however, some simulation parameters 
have to be different from those of the real world. The use of finite lattice spacing a and spatial 
extent L is unavoidable but its effects can be systematically reduced. It is assumed in the simu- 
lations reviewed in this article that up and down sea quarks are degenerate. The use of relatively 
heavy masses muci,sim for degenerate up and down quarks is much more problematic, because it 
could cause a large uncontrolled error by extrapolating lattice results to the physical mass m„j. It 
is therefore vital to simulate quark masses triud. sim, where chiral perturbation theory (ChPT) can be 
safely used as a guide for the chiral extrapolation. 

The staggered fermions are known to be computationally inexpensive [p, and led to a pre- 
cise determination of /k/ fn [§] from the MILC collaboration's simulations at '«Hd,sim^'«i/10 |^]. 
Their complicated flavor structure is however a serious obstacle to extensive calculations of more 
involved matrix elements, such as the K13 form factors. While simulations with other discretiza- 
tions were limited to relatively heavy quark masses, typically rnud.sim 

>ms/2, at the time of the 

last conference KAON 2005, recent algorithmic improvements now enable us to explore mud^saa 
comparable to that in the MILC's study. In any lattice studies, consistency between their data and 
ChPT is a crucial issue for a reliable chiral extrapolation. 

In this article, we first review recent progress on Jk/ fn in Sec. focusing on the latest update 
on the MILC's estimate and the status of studies with other discretizations. Section |3| is devoted 
to the Kit, form factors. We outline the calculation method and discuss the associated systematic 
errors. Finally, our conclusions are given in Sec. 0. 



As pointed out in Ref.|@], \Vus\ can be extracted from Ka and Tta decays through the ratio of 
their decay rates 



Radiative corrections parametrized by Ck.h and the muonic decay rates lead to an uncertainty of 
<0.2% in \Vus\- The determination of \Vud\ from super-allowed nuclear j8 decays is accurate at the 
impressive level of 0.05%. Therefore the main uncertainty in \Vus\ comes from the theoretical input 



Lattice QCD can, in principle, give a precise estimate of fx/ fn, since each decay constant 
is calculated from simple (and hence less noisy) two-point functions and uncertainties due to the 
lattice scale and renormalization are canceled in the ratio. The dominant error arises from the con- 



2. h/f: 




(2.1) 



hi fn- 



tinuum and chiral extrapolations. The original estimate \Vus\ =0.2236(30) in Ref.[Q] was obtained 
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with the MILC's estimate 7/^//;^ = 1.201 (8) (15) [Q] from their simulations using an improved stag- 
gered action (the so-called AsqTad action) at two lattice spacings a = 0.09 and 0.12 fm and with 
quark masses down to m„rf sim ~ m.s/10 

2.1 Update on result from staggered fermions 

The MILC collaboration has been steadily updating their simulations. One of the main im- 
provements in their latest report is that their simulations are extended to finer (a = 0.06 fm) and 
coarser lattice spacings (a = 0.1 5 fm) ||5|,^. 

It should be noted that simulations with the staggered quark action have the following theoret- 
ical and technical complications. By construction, a single staggered field describes four species 
of quark. This degree of freedom is called "taste". In simulations with single flavor (two degener- 
ate flavors) of quarks, gauge configurations are generated by taking the fourth (square) root of the 
quark determinant in the Boltzmann weight, e.g. 

Zn,=i= J[dU]det[D]'/^exp[-Sg], (2.2) 

where D is the Dirac operator for the four-taste staggered quark, Sg is the lattice gauge action 
of choice, and dU represents the path integral over the gauge fields. It is still actively debated 
whether the non-local Dirac operator corresponding to the rooted determinant leads to the correct 
continuum limit |^. In addition, the explicit taste symmetry breaking at finite lattice spacings 
makes calculations of matrix elements complicated. 

In the MILC's latest analysis [||, ^, they fit the quark mass and lattice spacing dependence of 
the pseudo-scalar meson masses and decay constants simultaneously using formulas from the so- 
called staggered ChPT where effects due to the taste symmetry breaking are taken into account. 
Analytic and chiral logarithmic terms at NLO and a part of analytic terms up to NNNLO are 
included into their fitting function to achieve a reasonable value of ;t^^/dof. Their two observations 
increase the reliability of their chiral and continuum extrapolations: 

• their fit curve in the continuum Umit exhibits a curvature towards the chiral Umit as expected 
from NLO ChPT; 

• they obtain low-energy constants (LECs) L4 = 0.1(4) and L5 = 2.0(4), which are consistent 
with a phenomenological estimate L4 = 0.0(8) and L5 =2.3(1) [^. 

They obtain their latest estimate 

^ = 1.208(2)(+7/-14), (2.3) 

where the first error is statistical and the second is systematic, and obtain /;;: = 128.6(0.4) (3.0) MeV 
and /a- = 155.3(0.4)(3.I) MeV, which are in good agreement with experiment. The statistical error 
is remarkably reduced from their previous estimate in Ref.[Q]. The uncertainty in fK/fn is now 
dominated by the systematics of the combined chiral and continuum extrapolation, which might 
be difficult to improve drastically without extending their simulations to a much wider range of 
fHud,sim and a. Independent calculations with different fermion discretizations are highly required 
to reduce the systematic uncertainties and to confirm that the rooted staggered theory has the correct 
continuum Umit. 
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2.2 Status of studies with otiier discretizations 



Unquenched simulations with other fermion discretizations also have a long history, leading 
up to recent studies with the following fermion actions: 

• (improved) Wilson fermions 

This traditional formulation is computationally cheap compared to chiral fermions (see be- 
low) and useful to simulate large and fine lattices. Its main drawback is the explicit chi- 
ral symmetry breaking at finite lattice spacing, which may distort the chiral behavior of 
the decay constants The clover action is an improved formulation by removing lead- 
ing 0{a) discretization errors. These discretizations are employed in recent simulations in 

Refs.[|r],|g,0]. 

• twisted mass Wilson fermions 



This is a variant of the Wilson fermions with the so-called twisted mass term [14], which sim- 



plifies the mixing pattern in the renormalization of lattice operators remarkably with compu- 
tational costs comparable to Wilson fermion simulations. This mass term, however, leads to 
the explicit breaking of parity and isospin symmetry. Its effects have to be studied carefully, 
as in large-scale simulations by the ETM collaboration [|T^]. 

chiral fermions 



With the five dimensional domain- wall formulation [16], chiral symmetry is restored in the 



limit of infinitely large size in the fifth dimension. It is however Ls/a times costly with 
respect to the above mentioned Wilson-type fermions. The RBC/UKQCD collaborations 
simulate three-flavor QCD with Ls/a = 16, which leads to the additive quark mass renor- 



malization of a few MeV []17|]. The (four dimensional) overlap fermions [ ]18| ] are even more 
computationally demanding. However, it has almost exact chiral symmetry and hence is use- 
ful for phenomenological applications where chiral symmetry plays an important role. The 
JLQCD collaboration has started large scale simulations in two-flavor QCD [p^]. 



The simulation cost for the above formu- 
lations with the commonly used Hybrid Monte 
Carlo (HMC) algorithm [EO] rapidly increases 



sqrt(2)M„ 



decreases: it scales as oc mj^-^ 
This is why previous simulations on relatively 
large and fine lattices |||, H, ||, || |6H7l 



KS; MILC, 2001- 

Wilson, SESAM/TzL, 1998 

clover; UKQCD, 2002 

clouer; CP-PACS, 2002 

clover; JLQCD, 2003 

clover; CP-PACS+JLQCD, 2003- 

domam-wall; RBC. 2005 

Wilsori; qq+q, 2003 

Wilsori; Del Debbio et al., 2007 



28| ] are limited to heavy quark masses typically 



wJud^sim^ as shown in Fig. |l|. However, re- 
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cent algorithmic improvements [ p9| , enable 
us to simulate much smaller values of m„j sim, 
which are now comparable to those in the MILC's 
simulation with the staggered fermions. 

In Fig. ^, we plot the pion decay constant 
obtained with Wilson-type and chiral fermions [15, 19, 26, 2^]. We observe a reasonable agreement 
among the data suggesting that discretization errors are not large in this plot. More importantly, 



Figure 1: Region of pion mass simulated in large- 
scale unquenched calculations in two-flavor (shaded 
band) and three-flavor QCD (solid band). 
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data at m„j sim ^ from recent simulations 
show a curvature toward the chiral limit as sug- 
gested by the chiral logarithm at NLO, whereas 
the curvature is not clear at heavier mud,sim- 

While 

data at small m„d sim are subject to effects due to 



finite lattice volume, the ETM collaboration JlSQ 
observe that their data with finite volume correc- 
tions [||] are described by the NLO ChPT for- 
mula reasonably well. They obtain 

F = 121.3(7) MeV, [4 =4.52(6), (2.4) 



which are consistent with lattice estimates of F 
from the MILC's simulation in /7-regime [0] and 
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Figure 2: Pion decay constant as a function of pion 
mass squared. Dashed line shows a linear fit to open 
symbols. 



JLQCD's one in £-regime [32], and with a phenomenological estimate of I4 [33]. This suggests 
that recent simulations with Wilson-type and chiral fermions are now exploring mud,sim sufficiently 
small to make contact with NLO ChPT. 

Recent estimates of /k / fn in three-flavor QCD 

are collected in Fig. |. The CP-PACS and JLQCD ^ 
collaborations obtain a slightly smaller result than ^ 



CP-PACS/JLQCD, 2006 
PACS-CS, 2007 
RBC/UKQCD, 2007 



1.05 1.10 1.15 1.20 1.25 1.30 



Others [p4[], probably because their simulations are 
limited to m„rf sim^'«i/2 and fji is overestimated 
due to the lack of the chiral logarithm. The PACS- 
CS collaboration employs the clover fermions with 

the Liischer's domain-decomposed HMC [pO|]. A ^ 

good agreement of their [0 and RBC/UKQCD's pigure 3: Recent unquenched estimate of fK/U 
estimates [|T^] with Eq. (2.3) is very encouraging, 

though their simulations are still on-going and/or the quoted error is statistical only. These are not 
enough mature to be used to derive an world average, and we simply quote Eq. ( |2.3| ) as the current 
best estimate of /k/ fn- It is, however, worth emphasizing that estimates of fx/ fn are expected to 
be improved remarkably in the near future by on-going simulations with the Wilson-type and chiral 
fermions by various groups. 



3. K12 form factor 

The Kit, decays provide a precise determination of | V„i | through 



T{K^Tilvi) 



192713 



|/+(0)|25EW(1+AEM + Asf;(2)), 



(3.1) 



where C is the Clebsh-Gordon coefficient equal to 1 ( 1 / \/2) for neutral (charged) kaon decays. 
The short- and long-distance radiative corrections, denoted by ^ew and Aem> and SIJ (2) breaking 
corrections are theoretical inputs, whereas the decay rate F and the phase space integral / 

are determined from experimental measurements. The uncertainties in IV,,.,! due to these inputs are 



well below 1% [36]. 
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Nf quark action a[fm] L[fm] Mps[MeV] 


JLQCD [grj 

RBC [0] 


2 clover 0.09 1.8 > 600 
2 domain-wall 0.12 1.9 > 490 


Fermilab/MILC 
RBC/UKQCD 1 


/H 
44 


PQCD 1 
] 


43 


] 


3 AsqTad (rf=clover) 0.12 2.5 > 500 
3 domain-wall 0.12 1.9,2.9 > 300 



Table 1: Simulation parameters in unquenched lattice calculations of /+(0). The clover fermions are used 
for valence down quarks in Ref . . 

The dominant uncertainty of the present estimate of \Vus\ therefore arises from theoretical de- 
termination of the normalization of the vector form factor /+ (0) defined from the A" — > 71 matrix 
element {n{p')\sYnu\K{p)) = {p + p')^U{q^) + {p- 
ing correction [|^] in the chiral expansion 



- p')^ /_ [q^], where = {p — p'Y. The lead- 



/+(0) = l+/2+/4 + O(/ 



(3.2) 



is practically free of uncertainties (/2 = —0.023), because any poorly known LECs do not appear 
in /2 thanks to the AdemoUo-Gatto theorem 

However, the higher order correction /4 contains LECs in the chiral Lagrangian both at 0{p'^) 
and 0{p^). A phenomenological estimate /4 = —0.016(8) based on the quark model was obtained 



by Leutwyler and Roos (LR) |39], and has been used in previous determinations of \Vus\- There 
has been remarkable progress in studies based on ChPT, where the evaluation of the tree-level 



, = 0.05 



contribution with LECs in the 0{p^) Lagrangian is the most crucial issue [4C]. Recent estimates 
ranging from /4 = — 0.007(9) to -1-0.007(12) are slightly larger than the LR estimate due to a (partial) 
cancellation between loop and tree-level contributions. 

Lattice QCD can provide a non-perturbative 
determination of /+(0), namely a«cf higher or- 
der contributions. Unquenched calculations per- 
formed so far are listed in Table |T[ These studies 
basically follow the strategy proposed in the first 
calculation in quenched QCD [p5[], which is out- 
lined below. 

The first step is to calculate the scalar form 
factor fo=U + {q^/{Ml-Ml))f^ from three 
point functions, e.g. 



{0,{t")\V^{t')\OK{t)^), (3.3) 




Figure 4: Double ratio Eq. (3^) as a function of t' 
with fixed f =4 and f" = 28. Data are from Ref. I 



where Oj^^f^^{t) is the interpolation operator for 
pion (kaon) and V^{t) is the vector current at the timeslice t. With sufficiently large temporal 
separations t" —t' and t'—t, C^^{t,t' ,t") is dominated by the ground state contribution, which is the 
matrix element (7r|V^|A') times unnecessary factors, such as the damping factor e^'^'^^' These 
factors are canceled in the so-called double ratio [fl6|]. For instance, a double ratio 



Cf'^{t,t',t")Cf{t,t',t") ^ {MK + M^f 
Cf''{t,t\t")Cl^{t,t',t") ^ AMkM^ 



l/o(^n 



(3.4) 
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KTeV 
ISTRA+ 
NA48 
KLOE 

FlaviaNet (fit) 



Becirevic et al.; Nj=0, clover 
JLQCD; N,=2, clover 
RBC: N|=2, domain-wall 
RBC/UKQCD; N,=3, domain-wall 



can be determined precisely as shown in Fig. ^, and it gives /o at q^^^ = {Mk—Mj^Y with an 
accuracy well below 1%. We can calculate /o at T^^max from different double ratios proposed in 
Refs. [45, 42], which however involve three-point functions with nonzero meson momenta and 
hence are much noisier than the ratio Eq. (3.4). 

Then, we interpolate /o to = 0. The q^ 
dependence is parametrized using the monopole 
ansatz /o(^^) =/o(0)/(l — Ao <7^) or polynomial 
forms up to quadratic order /o(<7^)=/o(0) + Ao^^ + 
Aq^^. These forms are also employed in analy- 
ses of experimental data. It turns out that the 
choice of the interpolation form does not cause 
a large uncertainty in the unquenched studies, 
since an accurate estimate of /o(<?max) avail- 
able near the interpolation point q^ = 0. It is also 
encouraging to observe in Fig. |5| that Xq from lat- 
tice studies shows a reasonable agreement with 
experimental measurements [^]. 

Finally, /+(0) =/o(0) is extrapolated to the 
physical quark masses m„d and m^. In all un- 
quenched calculations, a ratio motivated by the 
AdemoUo-Gatto theorem 




0.4 

Mps' [GeV'] 



Figure 5: Lattice estimates of slope A p [p3[ , pi , 42 
H] together with experimental values 



R 



/+(0)-l-/2 



(3.5) 



can be fitted to a rather simple polynomial form 



R = co + c,{Ml + Ml). 




RBC+UKQCD: N,=3, domain-wall 
1 + (ChPT) -1- NNLO aniy 
U NNLOanly 



(3.6) 



0.3 0.4 

2 2 

5 [GeVi 



0.5 



It is possible that, since most simulations are lim- 
ited to heavy quark masses fw„rf sim the 



Figure 6: Vector from factor as a function of 
'«!(i/,sim- The solid line shows a fit curve Eq. (3.6). 
NNLO (and higher order) chiral logs vary smoothly The NLO term /z is subtracted in the dashed line; 

. ^, . . J • ^ J 1. namely the difference between two lines shows /?. 

m this region and are well approximated by the 

analytic terms. 

In order to get an idea about how small rriud. sim is needed to see the chiral logs in /+ (0) clearly, 
data from the RBC/UKQCD's study is plotted a.s a. function of ifiud sim 

Fig. g The NLO chiral 

log /2 rapidly increases at m„rf sim ^ nis/l. This suggests that precise lattice data in this region are 
essential for a reliable chiral extrapolation compatible with the existence of the chiral logs. 

We note that the error of f+{0) may rapidly increase with decreasing niud, sim, because of longer 
auto-correlations of gauge configurations and larger q^^^ for the q^ interpolation of fo{q^). In fu- 
ture studies at small m,,^ sim> therefore, it is advisable to employ improved measurement methods, 
such as the all-to-all quark propagators to improve the accuracy of fo{q^) p8|]. The twisted bound- 
ary condition, which enables us to explore ~0 [49], and a model independent parametrization of 
the q^ dependence of /o [ 50 ] are useful to reduce systematic uncertainties due to the q^ interpola- 
tion. 
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Figure ^ shows recent lattice estimates of /+(0). The RBC/UKQCD collaboration confirms 
that finite volume corrections at L ~ 2 fm are small down to ifind ^lu^ Tunis /A. This observation is 
encouraging since other unquenched studies are conducted with similar or larger lattice sizes. The 
nice consistency among lattice results may suggest that discretization and quenching errors are not 
large. 

All lattice results are in good agreement with 
the LR value. We note that, however, estimates 
from ChPT are slightly higher due to the NNLO 
loop contributions. Therefore, the agreement be- 
tween lattice and the LR value has to be exam- 
ined carefully by precise lattice calculations at 



Leutwyler-Roos, 1984 



Bijnens-Talavera, 2003 
Jamin-Oller-Pich, 2004 
Cirigliano-Neufeld-PichI, 2004 
Cirigliano et al., 2005 



< 



ms/2, where chiral logs are expected 



to be seen clearly as discussed above. 



4. Conclusions 



N,=0 



N,=3 



Becirevic et al., 2005 
JLQCD, 2005 
RBC, 2006 

MILC, 2005 
RBC/UKQCD, 2007 



0.94 0.96 0.98 1.00 



f (0) 



Figure 7: Vector form factor /+(0) from phe- 
nomenological model (top panel), ChPT (middle 
panel) and lattice QCD (bottom panel). 



From the MILC's estimate of fx/ fn and the 
K^2/^^2 decay rates, we obtain \Vus\ =0.2226 
(+26/ — 15). The preliminary result of /+(0) = 

0.9609(51) from the RBC/UKQCD's study and |K.s/+(0)| = 0.21673(46) from the FlaviaNet 
working group [36] lead to \Vus \ =0.2255(13) which is consistent with the value quoted earlier. 
The latter estimate needs, however, further studies to increase the reliability of the chiral extrapo- 
lation of /+(0). 

For both fx/ fn and /+(0), we observe that precise data at sufficiently small ud quark masses, 
typically mud.sim^nis/2, are needed for a reliable chiral extrapolation. Thanks to recent algorithmic 
improvements, several groups have already started large-scale simulations in this region of rnuci,sim 
with different fermion discretizations. While their results are premature to be taken into account 
in the above estimates of \Vus\, lattice estimates of /k/ fn and Ki^, form factors are expected to be 



remarkably improved by these studies in the near future [51] 
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